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Abstract 

The idea of fractional exclusion statistics proposed by Haldane is applied to 
systems with internal degrees of freedom, and its thermodynamics is exam- 
ined. In case of one dimension, various bulk quantities calculated show that 
the critical behavior of such systems can be described by c = 1 conformal 
field theories and conformal weights are completely characterized by statis- 
tical interactions. It is also found that statistical interactions have intimate 
relationship with a topological order matrix in Chern-Simons theory for the 

fractional quantum Hall effect. 
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It is known that there are possibihties of existing particles which obey fractional statis- 
tics in lower dimensions. For example, in the fractional quantum Hall effect (FQHE) quasi- 
particles carry fractional charges, which may be described by Chern-Simons (CS) theory 
In this theory, as well as others, the dimension of space plays a crucial role because 
fractional statistics is defined on the basis of monodromic properties of wavef unctions. Re- 
cently, Haldane has proposed a new definition of fractional statistics independent of the 
dimension of space. It is based on the state-counting of many-body systems and can be 
considered as a generalization of the Pauli exclusion principle. Imagine that there is a A^- 
particle system with species ai,i = 1,2, ...,N. The A^-body wavefunction can be expanded 
by ith one-particle wavefunction by keeping the coordinates of other particles fixed. Let D^- 
be the dimension of such a basis. By observing that may be changed when a particle is 
added without changing the boundary and the size of the system, Haldane has introduced 
a statistical interaction Qa/s defined by 

where g^/s is a constant independent of {A^^}. The special cases of ga/3 = and ga/s = Sag 
correspond to free bosons and free fermions, respectively. Since Haldane's work, several 
authors have investigated it in detail 0-0. 

Wu [^] and Bernard and Wu ^ have studied the thermodynamics of single-component 
particles whose total number of states is given by 

]V.!(Bi-l)! • ^ ' 

with Di being a solution of eq.([l|), where the species a is regarded as the index i which 
labels the momentum ki. They have shown that gij is given by the two-body S matrix for 
Bethe ansatz solvable models. As simple examples in one dimension (ID), they calculated 
gij for the Bose gas with 5-function interaction and the Calogero- Sutherland model [Q. The 
former has non-zero off-diagonal elements of gij which they call mutual statistics, while the 
latter has only diagonal ones, gij = g6ij. Nayak and Wilczek P] have called such particles 



"(y'-ons" and examined low-temperature properties in arbitrary dimensions. It may be such 
g-ons that can be regarded as ideal particles obeying fractional statistics. 

In this Letter, we will generalize this approach to systems with internal degrees of free- 
dom. This has been briefly discussed in , and a step has been made forward in for the 
ID t/ — s> oo Hubbard model. First, we will obtain general formulae for thermodynamics, 
and then calculate bulk quantities in case of ID to find that critical behavior of such systems 
can be described by c = 1 conformal field theories (CFT). It will next be shown that there 
is intimate relationship between statistical interactions and a topological order matrix in CS 
theory for the FQHE ||. 

Let us consider a system of particles whose states are classified by the labels a and i, 
where a (= 1, 2, M) and i correspond to internal degrees of freedom and the momentum of 
particles, respectively. Denoting explicitly these two kinds of indices in eq.(|I|), the statistical 
interaction g^ipj for the multicomponent system reads. 



dNf,, 

which results in 



—QaiPj 5 (3) 



D^, = -Y^g^,p,Np,+Dl, (4) 
Pi 

where Nai {Dai) is the number of particles (holes) labeled by ai, and is the number 
of holes when there are no particles in the systems. Suppose that we have no particles in 
the system under consideration. Then there are only M kinds of free holes. Therefore, it 
is natural to choose D^^ = with = D^, which we will call a symmetric basis in this 
paper. Another important choice is = D^Sai, namely, a = 1 and a > 1 classify charge 
and internal degrees of freedom, respectively. This is referred to as a hierarchical basis. It 
is noteworthy that one can find the same bases also in a CS theory for the FQHE , as we 
shall see later. The distribution functions of particles and holes labeled by ai are defined as 



Charge for each component is also defined by ta = D^/ D^, which reduces to t = (1, 1, .., 1) 
and (1,0, ...,0) in the symmetric and the hierarchical bases, respectively. Note that the 
quantities (^) play a key role in the fractional exclusion statistics because they are general- 
izations of Fermi (or Bose) distribution function. In the following, we consider multicompo- 
nent systems whose total energy is given by the summation of bare single-particle energies, 
E = J2ai^ai^ai, whcre we set e°j = e^ta- It should be noted that many-body effects are in 
general incorporated in the distribution of Nai- 

Let us now examine the thermodynamics of a system whose number of states is given by 
eq.(0) with replacing i ai. Notice that eq.(^ is rewritten as 

It is now straightforward to investigate the thermodynamics following the method of Wu |^ 
and Bernard and Wu which is based on Yang and Yang's idea Thermal equilibrium 
is realized by minimizing the thermodynamic potential Q = E — TS — J2aPaNa, where 
Na = J2i and fi^ denote the total number of particles and the chemical potential for 
the species a, respectively. The entropy S is given hj S = \nW with W replaced i ^ ai in 
eq.(ll). The results are as follows. The free energy F at equilibrium is 

F = Y^ fiaNa -TJ2dI ln[l + exp(-e,i/T)], (7) 

" ai 

where the "dressed energy" e^j = T\n{Dai/Nai) = T \n{p^^^ / p^) is determined by 

eai-Tj2 9pjai ln[l + exp(-e/3j/T)] = 6°^ - (8) 
ft- 

where gai/sj = dai/sj — ^aipj- Though we here assume that there is no degeneracy for each 
state, we can also formulate the thermodynamics by introducing the corresponding weight 
factor in eqs.(|^ and (^) when each state has a certain degeneracy. The formulae (|^) and 
(^) describe thermodynamics of multicomponent systems obeying the fractional exclusion 
statistics (||). 
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Now consider low-temperature properties in the thermodynamic hmit, which is defined 
by Na oo and D° ^ oo proportionally 0. Hereafter, we restrict ourselves to ID systems 
for simplicity. Then eqs.(^~(^) read 

roo 

Pa\ka) =ta-Yl / dk/^gafBiK, kp)pp{kp), (9) 

p J —oo 

/oo 
dka\n[l+exp{-ea{K)/T)], (10) 

/oo 
dkfsgf3a{kf3, ka) ln[l + exp{-efs{kf^)/T)], (11) 

where p^^^ = + p'^^^. The total energy is given by e = E / = Y^a^a with Cq = 
IZodkaea{ka)pa{ka)- Assume gapi-ka,-kp) = ^0/3(^0,^/3) and el{-ka) = e°(/c„). Then at 
T = there exists a "generalized Fermi level" Qa which satisfies ea{ka) < for \ka\ < Qa 
and ea{ka) > for \ka\ > Qa- Hence Paika) = for |A;q,| > Qa, and Paika) for \ka\ < Qa 
and eaika) are given by the solutions of the following equations: 



Paika) =ta-2_^ / dkpQapika, kp)pp{kp), (12) 

dkpei3{kp)gfia{kp,ka). (13) 



Therefore, in the low-temperature limit, eq.(|TOD is expanded as 

F(T) = F(T = 0)--T2(27rD°)^— (14) 

6 Q, Va 

in terms of the Fermi velocity Va = ^'aiQa) I PaiQa)- The specific heat in low-temperatures 
in unit length is then given by C/{TL) = (tt/S) J^ai'^/'^a), where D° = L/{2-n) is used p. 
Hence, the critical behavior of ID systems whose statistical interaction is defined by eq.(^ 



can be described by M independent c = 1 CFT ||Tl[. The responses to "external fields" jl 
are also calculated as 

X./3 = 1^ = 2{ZV-'Z^)ap, (15) 

where Ua = Na/D^ = dkaPa{ka), Zai3 = ZapiQp) is a boundary-valued "dressed charge" 
Za/sikp) defined by [|l^ 



/Q-y 
dk^Za^{k^)g^p{k^,kii), (16) 

and a matrix V is defined by Va/3 = VaSa/3- From eq.(|l5D, the compressibility can be obtained 
as Xc = J2ai3taXai3ti3- The cxcitation spectrum in CFT hmit is classified as [jlO 

Ae = ^AniZ-^VZ-^An + AdZVZ^Ad (17) 

except for particle-hole excitations, where An (Ad) is a M-component vector of quantum 
numbers which change the number of particles (carry the large momentum). This expression 
also confirms that the critical behavior is described by M independent c = 1 gaussian CFT. 
Let us now consider a simple system with 

QaifSj = SijGa/S- (18) 

The choice (|l^) implies that the interaction in the momentum space is local, i.e., it takes 
a 5-function form. We refer to particles obeying (|T8|) as ideal multicomponent gases with 
fractional exclusion statistics which is characterized by the statistical interaction Gap- One 
can indeed find that eqs.(|^) and (|^) reduce to p^^K^) +1^/3 GapP/sik) = ta, which shows an 
asymmetry between particles and holes as discussed in refs. |lT2|,|5|,p!3[| . We will later discuss 



concrete models in ID which show fractional exclusion statistics (|T8|). From eq. (p^Sf) , the 
distribution of particles with species a at finite temperature is given by 

Paik)=J2G-'pik,T)tp, (19) 

where Gapik^T) = exp{ea{k) /T)6a/3 + Ga/3, and ea{k) is a solution of eq.(|TI|), i.e., 
n/3(e'''/^)*^'''»(l + e''5/^)'''3"-'^'5" = e^'^"'""^/^. We emphasize that the formula (|19D is a gen- 
eralization of the Fermi (or Bose) distribution function to a multicomponent ideal gas with 
the statistical interaction Gap- We here show simple examples in the symmetric basis. 
First, when Gap = g, distribution in eq.(|l9D becomes that of ideal fractional one-component 
particles, p = p^^\ which has been discussed in This formula further reduces to 

1 / (e*^^""^''/"^ — 1) for free bosons {g = 0) and to 1/ (e*^^""^^/-^ + 1) for free fermions {g = 1)- As 



another simple example of two-component system with Gu = G22 = 9 and 6*12 = G21 = f, 
we have p= {p^^'\ p^^'^) with g' = g + f. 

Consider now the T = case with general Gaf^. There appear M kinds of "generalized 
Fermi levels" Qi,Q2, ■■■,Qm- Define the ordering of a by T{a), such that Qr{i) > Qt(2) > 
... > Qr(M)- Then, in the region (5^(„) > \k\ > Qr{n+i), we find pr{a){k) = Y.p=i ^n(G')a/^^r(/3) 
for a = l,...,n and Pr(a){k) = for a = n + 1,...,M, where t„(G) is a n x n matrix 
whose elements are given by Tn{G)ai3 = Gr{a)T{f3) with a,P = l,2,...,n. In the absence of 
external fields one finds, by putting pa = pta, that the ground state is obtained by setting 
Qi = Q2 = ■■■ = Qm = Q- The dressed energy ea{k), distribution function pa and the 
density of particles ric = N^/ D^, where N^. = J2a ^ata, are calculated as 

ea{k) = j:G-'p{el{k)-pp), (20) 

Pa=T.G-'ptp^ (21) 

n, = 2Qiy, = J^taG-^tp. (22) 

These equations show that Fermi velocities are calculated as Va = Q = nc/(2h'), which are 
independent of their components. If we take e^{k) = as the bare spectrum, we then find 
the ground state energy eg = nj?/(24i/^) and the compressibility Xc = ^i''^ /ric- Hence, these 
bulk quantities are determined solely by v which does not involve full information of the 
matrix G. However, if we observe the excitation spectrum Ae = {v/ 4)^nG^n + v^dG~^A.d^ 
there shows up every element of the matrix G, owing to the fact that M kinds of Fermi 
velocities are the same (See eq.([T7D). Here we assume G^ = G. Note that this formula for 
the excitation is of a desirable form required by CFT, and we thus find conformal dimensions 



determined completely by the matrix G following the idea of finite size-scaling |]TT|. 

A remarkable point is that the above matrix G of statistical interactions has close re- 
lationship to a topological order matrix in the FQHE. To see this explicitly, it is more 
convenient to take e^{k) = k as the bare spectrum, which only contains right-going parti- 
cles. It is known that such chiral particles in ID may describe essential properties for edge 
states of the FQHE 0. In this case, the excitation spectrum of the present system has 



the form Ae = (f/4)AnGAn, which may be described by the holomorphic piece of CFT. 
Recall here that a certain hierarchy of the FQHE can be described by a multicomponent CS 
Lagrangian L = t^^^aa^Kapdi^aj^x in which the construction of a hierarchy is characterized 
by the topological order matrix K p|. An important point is that the excitation spectrum 



obtained here for the chiral model coincides with that for edge states [Tj] characterized by 



the topological order matrix K ]T5|]. Therefore we arrive at an instructive consequence that 
the matrix G of statistical interaction can be identified with a topological order matrix K 
in CS theory for the FQHE. 



For example, the topological-order matrix K used in ref. is written in the symmetric 
basis, 

M 



^=E(P"-2 + <5ai)CF + / (23) 



— ^ -r Ual)^ 



with an odd integer pi and even integers pa (a > 1), where / is a unit matrix and Ci,j is 
a M X M matrix defined by (C^/'')jj = 1 for i,j = a,a + 1, ...,M and = for otherwise. 
Substituting eq.(]23|) into eq. (^2]) , we have 



V 



1 

Pi Y 



(24) 



P2 

Pm 

Note that this quantity is nothing but the filling factor for the FQHE described by (p3|). 
It is known that the symmetric basis is converted into the hierarchical basis by the 
matrix W defined by W^p = Sa/3 — Sa+i,/3 as follows; t tW and K W^KW. In this 
basis, distribution functions (pT)) are easily calculated at zero temperature; pa = n^=i ^p, 
where = I]f=i G^li ^ with (3 x (3 matrix G'^^^ defined by G^-^'' = GM-i3+i,M-l3+j- 
particular, when pi = 2m + 1 and Pq, = 2 for a = 2,3, ...,M, we have the filling factor 
u = M/{2mM + 1) and the corresponding distribution functions p = {u, (M — 1)/M, (M — 
2)/{M — 1), 2/3, 1/2). This family corresponds to a fundamental hierarchy which has 
been observed in the FQHE 

Finally some comments are in order for concrete models in ID, which exhibit essential 
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properties of fractional exclusion statistics. It has been claimed that ideal particles with 
fractional exclusion statistics are realized by ID quantum models with interaction 
r7| , p!8[] and their multicomponent versions |T^. In fact, this fractional property has been 



shown explicitly for the Haldane-Shastry spin chain |^ and also for the Calogero-Sutherland 
model |T^p!3| , PH . For cases with internal degrees of freedom discussed here, we find that 
eqs.(|])~(|D indeed reduce to the Bethe ansatz equations for multicomponent Calogero- 
Sutherland model with a special coupling if we convert the expressions to those of the 



hierarchical basis. Hence, an ideal multicomponent gas with statistical interaction G may 
also be realized by multicomponent models. 

On the other hand, for interacting ID particles with general type of interaction, exci- 
tations cannot be regarded as ideal particles. Nevertheless, if we restrict ourselves to the 
CFT (low-energy) limit of such systems, the excitation can be well approximated by ideal 
particles with fractional exclusion statistics. For example, consider the ID Hubbard model 
in the CFT limit. There are two elementary modes called spinon and holon, which have 
different velocities Vg ^ Vh in general. Hence, it seems difficult to deduce the matrix G 
from eq.(^). However, these velocities are not universal, and therefore by choosing Vg = Vh 
hypothetically, we can conclude that statistical interaction G is given by G = I + aCP with 
a = {l/Kp — l)/2, where Kp (1/2 < Kp < 1) is the critical exponent (divided by 4) for 
the Akp oscillating part of the density-density correlation function This kind of inter- 
pretation for the elementary excitation is valid in general for the class of multicomponent 
Luttinger liquids. 

In summary, we have applied the idea of Haldane's fractional statistics to systems with 
internal degrees of freedom. The critical behavior of such ID systems can be described 
by multiple c = 1 CFT with conformal weights given by statistical interaction G. It has 
also been shown that there is intimate relationship between statistical interactions and a 
topological order matrix in CS theory for the FQHE. It is quite interesting and still open 
to study how the fractional exclusion statistics control the physics of correlated particles in 
higher dimensions. 
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